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Abstract 

We characterize the angular polyspectra, of arbitrary order, associated with isotropic fields defined 

on the sphere 5*^ = {{x, y, z) : x^ + y^ + z^ = l}. Our techniques rely heavily on group representation 

theory, and specifically on the properties of Wigner matrices and Clebsch-Gordan coefficients. The 

P^ , findings of the present paper constitute a basis upon which one can build formal procedures for the 

Qi^ ' statistical analysis and the probabilistic modelization of the Cosmic Microwave Background radiation, 

^ I which is currently a crucial topic of investigation in cosmology. We also outline an application to 

random data compression and "simulation" of Clebsch-Gordan coefficients. 
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^ 1 Introduction 
l^ ■ 

QQ The connection between probability theory and group representation theory has led to a long tradition 
\^ ' of fruitful interactions. A well-known reference is provided by |T^; see e.g. [H Section 40-41], [TT], [IB] . 
^D ' [Sj, [IS], [55], [31], and the references therein, for other relevant contributions. In this paper we shall 
r — - focus in particular on the profound connection between the probabilistic notion of isotropy, i.e. invariance 
^^ in law under the action of a group, and the representation theory of the group itself. One instance of 
OO , this connection is well-known, i.e. the celebrated Peter- Weyl Theorem, which allows the construction of 
spectral representations for isotropic random fields on homogeneous spaces of general compact groups, 
see [24. for a general construction and [53] , [55] for examples related, respectively, to the torus and the 
sphere. Our aim here is to use these representations in order to characterize random fields by means of 
a higher order spectral theory; in particular, one of our main goals will be to establish the link between 
C^ the so-called polyspectra (or higher order spectra) and alternative (tensor product and direct sum) 
representations of the underlying isotropy group. In particular, we shall provide a general expression for 
higher order spectra of isotropic spherical random fields in terms of convolutions of Clebsch-Gordan or 
Wigner coefficients. The latter where introduced in Mathematics in the XIX century for the analysis of 
Algebraic Invariants; they have since then played a crucial role in the development of Quantum Physics in 
the XX century (see for instance [31] for a comprehensive reference); their role in Group Representation 
theory will be discussed below, while more details can be found for instance in [32] . 

Our analysis may have an intrinsic mathematical interest, but it is also strongly motivated by applica- 
tions to Physics and Cosmology. Concerning the latter, the analysis of higher order spectra for isotropic 
spherical random fields is currently at the core of several research efforts which are related to the analysis 
of Cosmic Microwave Background (CMB) radiation data, see for instance [l5] for a general introduction 
and [TS] [ITl [SO] [21] for some references on the bi- and trispectrum. A general characterization of the 
theoretical properties of higher order angular power spectra can yield several insights into the statistical 
analysis of the massive datasets that are or will be made available by satellite experiments such as WMAP 
or Planck. For instance, the current understanding of the behaviour of the bispectrum for some simple 
physical models has already led to many applications ([H], [34) . [35j). aiming at obtaining constraints on 
nonlinearity parameters of utmost physical significance; needless to say, a proper understanding of higher 
order spectra can lead to more efficient statistical procedures and better constraints, which may help to 



solve some of the important scientific issues at stake in CMB analysis (primarily a proper understanding 
of the Big Bang inflationary dynamics, which is tightly linked with the CMB nonlinear structure, see 

[H], [5], i, uni). 

The relevance of the current results need not be limited to cosmological applications. Indeed, the 
analysis of spherical random fields has currently led to remarkable developments in the Geophysical and 
Planetary Sciences, and even in Medical Imaging, see for instance ([9], [27], [33]). Moreover, we shall 
show below how the relationship which we establish leads very naturally to some numerical algorithms for 
the estimation of Clebsch-Gordan and Wigner coefficients. The latter represent probability amplitudes 
of quantum interactions and as such a rich literature in Mathematical Physics has been concerned with 
recipes for their numerical estimation: our procedure lends itself to easy implementation and can be 
simply extended to very general compact groups, although in this paper we focus solely on 5*0(3). 

The plan of this paper is as follows: in Section [2] we introduce our general probabilistic setting and 
provide some preliminary notation and background material. In Section [3] we present some background 
material on representation theory, while in Section 2] and Section [5] we obtain our main results, including 
the explicit characterization of polyspectra. These results are applied in Section [S] to derive explicit 
expressions in some important cases (such as x^ random fields). Section [7] is devoted to further issues that 
we see as the seed for future research: they concern, in particular, the connection with the representation 
theory for the symmetric group and the Monte Carlo estimation of Clebsch-Gordan coefficients. 

In the subsequent sections, every random element is defined on an appropriate probability space 

2 General setting 

In this paper, we focus on real-valued, centered, square-integrable and isotropic random fields on the 
sphere 5^ = {(a;,y, 2) G M^ : a;^ + y^ + z^ = l}. A centered and square integrable random field T on S*^ 
is just a collection of random variables of the type T — {T{x) : x G S^} such that, for every x S 5"^, 
ET{x) = and ET'^{x) < 00. In the following, whenever we write that T is a field on 5*^, we will implicitly 
assume that T is real-valued, centered and square-integrable. From now on, we shall distinguish between 
two notions of isotropy, which we name strong isotropy and weak isotropy of order n (n > 2). 

Strong isotropy — The field T is said to be strongly isotropic if, for every fc G N, every xi, ...,Xk G S"^ 
and every g G 50(3) (the group of rotations in R^) we have 

{T{x,),...T{xk)}^{T{gxi),...T{gxk)} , (2.1) 

where = denotes equality in distribution. 

Weak isotropy — The field T is said to be n- weakly isotropic [n > 2) ii E\T{x)\^ < cxd for every x G 5^, 
and if, for every xi, ...,x„ G S"^ and every g G 5*0(3), 

E [T{xi) X • • • X T{xn)] - E [T{gxi) x • • • x T{gxn)] ■ 



The following statement, whose proof is elementary, indicates some relations between the two notions 
of isotropy described above. 

Proposition 1 1. A strongly isotropic field with finite moments of some order n > 2 is also n-weakly 
isotropic. 

2. Suppose that the field T is n-weakly isotropic for every n > 2 (in particular, £'|r(a;)|" < 00 for 
every n > 2 and every x G S'^ ) and that, for every k > 1 and every {xi, ...,Xk), the law of the vector 
{T (xi) , ...,T(xk)} is determined by its moments. Then, T is also strongly isotropic. 



Now suppose that T is a strongly isotropic field, and denote by dx the Lebesgue measure on S . Since 
the variance ET (x) is finite and independent of x (by isotropy), one deduces immediately that 
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T (x) dx 



< oo, 



from which one infers that the random path x —> T (x) is a.s. square integrable with respect to the 
Lebesgue measure. Then, it is a standard result that the following spectral representation holds: 



oo I -. 

T{x) = '^ X] "'™^'™ (^) ' where aim ^ T (x) Yim ( 

J n ™, 1 ^ S 



x)dx, 



(2.2) 



1=0 m 



and where the complex-valued functions {Yi^ : / > 0, m ~ — /, ..., ^} are the so-called spherical harmonics, 
to be defined below. The spectral representation (|2.2p must be understood in the L^{il x S"^) sense, i.e. 
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L — *oo 



T — 2_^ 2-^ '^ImXlr 



1=0 m=-l 



-0, 



L^iS"^) 



where L^{S^) is the complex Hilbert space of functions on S"^, which are square-integrable with respect to 
dx. If moreover the trajectories of T(x) are a.s. continuous, then the representation (|2.2p holds pointwise, 
i.e. 



lim {T{x)- 
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;=0 m=-l 



o-irnYim (x) > = for all X € S"^, a.s.-P, 



see for instance pQ or [3S]. The spherical harmonics {i^/m}„^_; ; are the eigenfunctions of the Laplace- 
Beltrami operator on the sphere, denoted by A52, satisfying the relation Ag2Yim = —l{l + 1)^^;™- These 
functions can be represented by means of spherical coordinates x — (0, (p) as follows: 



l2l + l{l-m) 



4tt (l + m)l 



Pirnicos 9) exp{imip) , for m>0 



Yim{0, (^) = (-l)"iV^(0, (^) , for m < , 0<e<TT, < ^ < 2tt 



where P/m(cos6') denotes the associated Legendre polynomial of degree I, m, i.e. 



Pimix) - (-l)"(l-a;2r/2 
m ^ 0,1,2, 



d" 



dx 
I , I = 0,1,2,3, 



Piix) , Pi{x) 
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2H\ dx^ 



{x'-iy, 



The random spherical harmonics coefficients {aim} appearing in (j2.2p form a triangular array of zero- mean 
and square-integrable random variables, which are complex-valued for m ^ and such that Eaim'WmJ = 
5\'5™'Ci, the bar denoting complex conjugation. Here, and for the rest of the paper, the symbol 5^ is 
equal to one if a = 6 and zero otherwise. We also write Ci = E\aim\ , ^ > 0, to indicate the angular 
power spectrum of T (we stress that the quantity Ci does not depend on m - see e.g. [J] for a proof 
of this fact). Observe that, by definition of the spherical harmonics, aim = (— l)'"ai_m. Note also that 
a convenient route to derive (|2.2p is by means of an appropriate version of the stochastic Peter-Weyl 
theorem - see for instance [5] or [24], as well as Section [3T] below. 

Observe that the representation (|2.2p still holds for fields {T{x)} that are not necessarily isotropic, but 
such that the random path x — > T{x) is P-a.s. square integrable with respect to the Lebesgue measure 
dx. Indeed, if the last property holds, then one has that, P-almost surely. 



lim / T{x) - ^ ^ aimYim (x) 

^°^JS^ \ l=Om=-l J 



dx^O. 



(2.3) 



In this case, however, none of the previously stated properties on the array {«;,«} holds in general. By 
an argument similar to those displayed above, a sufficient condition to have that x —> T{x) is P-a.s. 
Lebesgue-square integrable is that sup^g52 ET{xY < oo. 

The next result, that we record for future reference, is proved in [3]. 

Proposition 2 Let T be a centered, square-integrable and strongly isotropic random field. Let the coef- 
ficients {aim} be defined according to i2.2\) . Then, for every l,m, one has that E\aim\ < oo. Moreover, 
for every I > 1, the coefficients {o/o, •■•, a;;} are independent if, and only if, they are Gaussian. If the 
vector {aiQ, ...,au} is Gaussian, one also has that ^(aim) and 'i (aim) are independent and identically 
distributed for every fixed m — l,...,l (^{z) and 3 (z) stand, respectively, for the real and imaginary 
parts of z). 

The following result formalizes the fact that, in general, one cannot deduce strong isotropy from weak 
isotropy. The proof makes use of Proposition [TJ 

Proposition 3 For every n > 2, there exists a n-weakly isotropic field T such that T is not strongly 
isotropic. 

Proof. Fix / > 1, and consider a vector 

'^ni 7 ^m — t , . . . , t , 

of centered complex-valued random variables such that: (i) feg is real, (ii) b-m = (~1) bm {m — 1, ..., I), 
(iii) the vector {6o, ..., bi\ is not Gaussian and is composed of independent random variables, (iv) for every 
k — 1, ...,n, the (possibly mixed) moments of order k of the variables {b^, ...,6;} coincide with those of a 
vector {flo, ..., a;} of independent, centered and complex- valued Gaussian random variables with common 
variance Ci and such that oq is real and, for every m = l,...,l, the real and imaginary parts of «„ are 
independent and identically distributed (the existence of a vector such as {6o, ■••,&;} is easily proved). 
Now define the two fields 

I I 

T (x) ^ 51 ^rriYlm {x) and T* (x) = ^ OrnYlm (x) . 
m— — l m^—l 

By Proposition [21 T* is strongly isotropic, and also n-weakly isotropic by Proposition[TJ By construction, 
one also has that T is n-weakly isotropic. However, T cannot be strongly isotropic, since this would 
violate Proposition [2] (indeed, if T was isotropic, one would have an example of an isotropic field whose 
harmonic coefficients {bo, ...,6;} are independent and non-Gaussian). ■ 

In what follows, we use the symbol A ^ B to indicate the Kronecker product between two matrices 
A and B. Given n > 2, we denote by II (n) the class of partitions of the set {1, ..., n}. Given an element 
TT S n (n), we write tt = {6i, ...,5fc} to indicate that the sets bj C {1, ...,n}, j = 1, ..., k, are the blocks of 
TT. The blocks of a partition are always listed according to the lexicographic order, that is: the block bi 
always contains 1, the block 62 contains the least element of {1, ..., n} not contained in 61, and so on. Also 
the elements within each block bj are written in increasing order. For instance, if a partition vr of {1, ...,5} 
is composed of the blocks {1, 3} , {5, 4} and {2}, we will write tt in the form tt ~ {{1, 3} , {2} , {4, 5}} . 

Definition A. (Al) Let the field T admit the representation (|2.2p . and suppose that, for some n > 2, 
one has that E |a;,„|" < cxd for every /, m. Then, T is said to have finite spectral moments of order n. 
(A2) Suppose that T has finite spectral moments of order n > 2, and, for I > 0, use the notation 

ai. = iai_i,...,aiQ,...,au) . (2.4) 

The polyspectrum of order n — 1, associated with T, is given by the collection of vectors 

Si,. ..I,, = E [ai,, ® ai^. (g) • • • «) «;„.] , (2.5) 



where < li,l2, ■■■Jn- Note that the vector S'/j...;„ appearing in (|2.5p has dimension (2Zi + 1) x • • • x 
(2/„ + l). 

(A3) Suppose that T has finite spectral moments of order n >2. The (mixed) cumulant polyspectrum 
of order n — 1, associated with T, is given by the vectors 



7r={fci,...,bfc}en(n) 

where < Zi, l2, ■■■, In, and, for every block bj = {ii, ..., Zp}, we use the notation 



(2.6) 



E[(^,^b,ai^] =E 



ai- ® ■ ■ ■ ® ai- 



(recall that we always list the elements of bj in such a way that ii < ■ ■ ■ < ip) ■ Plainly, the vector Sf ; 
in (|2.6p has also dimension {2li + 1) x • • • x (2Z„ + 1) . 

Remark. Suppose that T has finite spectral moments of order n > 2. Then, by selecting frequencies 
h ^ I2 ^ ■ ■ ■ — h = I ^ Qj one obtains that 
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^'l...l ■■-SlAn)= J2 i-ir'{k-iy.E (a,) 

"i^, 7i-={fci,...,6fc}en(n) 



(!S>\bi\ 
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(2.7) 



where \bj\ stands for the size of the block bj, and we use the notation 



(ai.) 



0\bj\ 



ai (g) • • • (g) a; 

^ V ' 

I bj I times 



3 Preliminary material 

3.1 Representation Theory for 5*0 (3) 

We start by reviewing briefly some background material on the special group of rotations SO{3), i.e. the 
space of 3 X 3 real matrices A such that A' A = I3 (the three-dimensional identity matrix) and det(y4) = 1. 
We first recall that each element g G 50(3) can be parametrized by the set (</?, 1?, iIj) of the so-called Euler 
angles {0 < f < 27r, 0<??<7r, O<'0< Stt); indeed each rotation in M.^ can be realized sequentially as 



A - Aig) = Ri^J,^,^) = R,i^)R,{i9)R,i^) 



(3.8) 



where Rz{ip),Rx{^), Rz{ip) G S0{3) can be expressed by means of the following general definitions, valid 
for every angle a. 
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cos a 


— sin a 
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sin a 


cos a 





, i?.(a) = 





cos a 


— sin a 








1 / 
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sin a 


cos a 



The representation p.8p is unique except for i9 = or z? = tt, in which case only the sum ip + tp is 
determined. In words, the rotation is realized by rotating first by ip around the axis z, then rotating 
around the new x axis by ??, then rotating by ip around the new z axis. It is clear that the first two 
rotations identify one point on the sphere, so the whole operation could be also interpreted as moving 
the North Pole to a new orientation in S^ and then rotating by ip the tangent plane at the new location. 
In these coordinates, a complete set of irreducible matrix representations for 5*0(3) is provided by 
the Wigner's D matrices D\ip,d,ip) = {^D^^^{ip,d,tp)]^^^_^ ^ of dimensions {21 -|- 1) x {21 + 1) for 
I = 0, 1,2, ...; we refer to classical textbooks, such as [32], [2] or [10], for any unexplained definition or 



result concerning group representation theory. An analytic expression for the elements of Wigner's D 
matrices is provided by 

Dl^i^P, d, if) = e-*"'^^^ We^™^ m,n = -{2l + l), ..., 21 + 1 

where the indices r?i, n indicate, respectively, columns and rows, and 



X 



E(-i)'^ 



k 



k\(l — m ^ ky.{l — n — ky.(m + n + A:)! 



and the sum runs over all k such that the factorials are non-negative; see [3T1 Chapter 4] for a huge 
collection of alternative expressions. Here we simply recall that the elements of D\il;,d,ip) are related 
to the spherical harmonics by the relationship 



Att / Att 

DUv,^,^') = (-i)'"y^^^r,-™(^,^) = y ^7^^r™(^,^) ■ (3.9) 

In other words, the spherical harmonics correspond (up to a constant) to the elements of the "central" 
column in the Wigner's D matrix. Such matrices operate irreducibly and equivalently on (2Z -I- 1) spaces 
(the so-called isotypical spaces) , each of them spanned by a different column n of the matrix representation 
itself. The elements of column n correspond to the so-called spin n spherical harmonics, which enjoy a 
great importance in particle physics and in harmonic expansions for tensor valued random fields. In this 
paper, we restrict our attention only to the usual n — spherical harmonics, which correspond to usual 
scalar functions. 

Remark. By exploiting relation (|3.9p , it is not difficult to show that the usual spectral representation for 
random fields on the sphere, as given in (|2.2[) . is just the stochastic Peter- Weyl Theorem on the quotient 
space S^ = SO{3)/SO{2). Indeed, by the stochastic Peter-Weyl Theorem (see e.g. [24]) we obtain, for 
any square integrable, isotropic random field {T(g) : g G 50(3)} 



T{g) = T{ip, ^,^)=Y.Y. ^IrnnJ^T^Dl^i^, d, ^) 



. 87r2 



where dg is the Haar (uniform) measure on 5*0 (3) with total mass Stt^. Now if we consider the restriction 
of T(5) to S^ = S'0(3)/5'0(2), denoted by Ts2{ip,d), we deduce that 



airnn = / Ts2{g)\ —-^D^„{g)dg 

Jso(3) V 8tt^ 



SO(3) 

2ir 






Ts.{ip,d){ I e^'^dVh/^^rfLWe-^'^sinMy.d^ 



2/ -Kl , 



^ T5.(^, z?)<5,",(27r)y ^-^dL„(i?)e-™'^ sinM^di? , 

the second equality following from the fact that T52 [g] is constant with respect to i}). We can thus conclude 

that 

f for n 7^ 

"" [ ^/2^ialm for n — ' 
where the array {aim} is defined by 



3.2 The Clebsch-Gordan matrices 

It follows from standard representation theory that we can exploit the family |Z3'} to build 

alternative (reducible) representations, either by taking the tensor product family |Z3'i (g) -D'^}, , or by 

considering direct sums < ©;Lu^-; 1^ i ■ These representations have dimensions 

(2;i + l)(2/2 + l)x {2h + l){2l2 + 1) 
and are unitarily equivalent, whence there exists a unitary matrix Ci-^i^ such that 

{D^^^D'^} = G,,, {®tt-i,\D^}c;^h ■ (3-10) 

The matrix C/^/j is a {{2li + 1)(2^2 + 1) x (2Zi + 1)(2^2 + 1)} block matrix, whose blocks, of dimensions 
(2^2 + 1) X {21 + 1), are customarily denoted by C^ , w , mi = —li, ..., li; the elements of such a block 
are indexed by m,2 (over rows) and m (over columns; note that m = — (2^ + 1), ..., 2^ + 1). More precisely. 



C, 



hh 






(3.11) 

%l= — ll,...,li;l = \l2—li\:...,l-2+ll 
h{mi)l2. ~ I hmihm2)m2 = -h,---,l2;m=-l,...,l ' \'^-^^) 

Remark. The fact that the two matrices D'^ D^^ and ®,^t/ , ,D^ have the same dimension follows 
from the elementary relation (valid for any integers ^1,^2 > 0): 

J2 (2/ + l) = (2^1 + 1) (2^2 + 1). (3.13) 

l=\l2-h\ 

By induction, one also obtains that, for every n > 3, 

E E •■• E (2A„_i + l) = n(2^, + l), (3.14) 

Ai = |b-/i| A2 = |;3-Ai| A„_i = |/„-A„_2| J = l 

for any integers ^1, ..., Z„ > (relation (|3.14p is needed in Section [5l2|) . 

The Clebsch-Gordan coefEcients for 5*0(3) are then defined as the collection {C'™,^ / m } of the the 
elements of the unitary matrices Ci-^i^. These coefficients were introduced in Mathematics in the XIX 
century, as motivated by the analysis of invariants in Algebraic Geometry; in the 20th century, they have 
gained an enormous importance in the quantum theory of angular momentum, where Cj"^ ; ^ repre- 
sents the probability amplitude that two particles with total angular momentum li,l2 and momentum 
projection on the z-axis ttii and TO2 are coupled to form a system with total angular momentum / and 
projection m (see e.g. |18|). Their use in the analysis of isotropic random fields is much more recent, see 
for instance [15j and the references therein. 

Remark {More on the structure of the Clebsch-Gordan matrices). To ease the reading of the subsequent 
discussion, we provide an alternative way of building a Clebsch-Gordan matrix C;^/^, starting from any 
enumeration of its entries. Fix integers ^i, ^2 ^ such that ^i < I2 (this is just for notational convenience), 
and consider the Clebsch-Gordan coefficients {Cji"^ j ^ } given in (|3.1ip - (|3.12p . According to the above 
discussion, we know that: (i) —li < mi < k for i = 1, 2, (ii) h — li<l<li +I2, (hi) —l<m< I, and (iv) 
the symbols {li,mi,l2,m,2) label rows, whereas the pairs {l,rn) are attached to columns. Now introduce 
the total order <c on the "column pairs" {l,m), by setting that {l,m) <c (^',"i'), whenever either / < V 
OT I = V and m < m' . Analogously, introduce a total order -<r over the "row symbols" (Zi, toi, I2, rn2), by 



setting that (?i,tii, Z2, '712) <r (^ii "z'j^, Zj, '^12), if either mi < m'^^ or mi = m'^ and m2 < m'2 (recall that 
li and I2 are fixed). One can check that the set of column pairs (resp. row symbols) can now be written 
as a saturated chaiiij with respect to -<c (resp. -<r) with a least element given by {I2 — li,~ {h ~ h)) 
(resp. {li, —li, ^2, "^2)) and a maximal element given by {I2 + ^i, h + ^1) (resp. (Zi, ^i, Z2, ^2))- Then, (A) 
dispose the columns from west to east, increasingly according to ^c, (B) dispose the rows from north to 
south, increasingly according to ^r- For instance, by setting li = and I2 > 1, one obtains that C'l-^i^ 

is the {2I2 + 1) X (2^2 + 1) square matrix I Cqq™^^^ > with column indices to — —{2I2 + 1), ..., (2^2 + 1) 

and row indices TO2 = —(2^2 + 1), ■•■, (2^2 + 1) (from the subsequent discussion, one also deduces that, in 
general, CgJJJ^^^ = 5^5^^). By selecting li = I2 = 1, one sees that Cn is the 9x9 matrix with elements 
Cimiima (^^^ TOi, TO2 = —1,0, 1; Z = 0, 1, 2, TO = —I, ..., I) arranged as follows: 



/ ^0,0 ^1,-1 



-(0,0 
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^0,0 
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'-'I, -1:1,1 



C2, — 2 ^2.-1 

1.-1:1,-1 '-'1,-1: 
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/--2,-l 
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1,0;1,-1 
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'-'1,1 
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p2, 
'-'1,1 
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'-^1,1: 



'-'1,-1 
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1,1 '-'1,1;1,1 / 



l.-l\ 



Explicit expressions for the Clebsch-Gordan coefficients of 5*0(3) are known, but they are in general 
hardly manageable. We have for instance (see [3T], expression 8.2.1.5) 



Iimil2m2 



(-l)'i+'^'+"V2^3 + l 



(;i+;2-?3)!(^l-^2-f;3)!(^l"^2 + /,3)! 



(^3 + msYHs ~ m^) 



{h + h + h + iy. 

1/2 



E 



(Zi + TOi)!(/i - TOi)!(/2 -I- TO2)!(Z2 - ™2)!_ 

{-lYih + I3 + mi - zy.ih - TOi + z)! 



1/2 



Zl{h + h-h- Z)\{h + TO3 - Z)\{ll -h- TO3 -f z)\ ' 



where the summation runs over all z's such that the factorials are non-negative. This expression becomes 
much neater for to,i = 7712 = m^ = 0, where we have 



, for ?i -I- /2 -I- '3 odd 



W30 

'-'ZiOZ20 



(-1) 



'l + '2-'3 



V2h + l[(h+l2+l3)/-2] 



[{ll+l2-l3)/2]\[(ll-l2+l3)/2]\[(-l 

for I1+I2 + I3 even 



l+/2+/3)/2]! \ 



+i2-hy.{ii-i2+hy.{-h+i2+hy. 

(h+h+h + iy 



1/2 



The coefficients, moreover, enjoy a nice set of symmetry and orthogonality properties, playing a crucial 
role in our results to follow. From unitary equivalence we have the two relations: 



^I\'n\\l2m2 I\ra\l2m2 / m ' 

Tii,m2 

'^Ixmxl2m2^lim'^l2ra'^ ~ ''miOm2 , 



(3.15) 
(3.16) 



l,m 



^ Given a finite set A = {aj : j = 1, ..., N} and an order -< on A, one says that A is a saturated chain with respect to -< 
if there exists a permutation tt of {!,..., N} such that 



«7r(l) -< «7r(2) -< ■■ ■ -(a. 



ir(iV-l) -< 0.tt(N) 



In this case, a^(i) and a^rpf\ are called, respectively, the least and the maximal elements of the chain (see 1281 p. 99]) 



in particular, (j3.15p is a consequence of the orthogonality of row vectors, whereas (I3.16|) comes from the 
orthogonality of columns. Other properties are better expressed in terms of the Wigner's coefficients, 
which are related to the Clebsch-Gordan coefficients by the identities (see [21], Chapter 8) 

^ >- ^> /2j—ri h-rnil2~m2 1«>--L'j 



h 


h h \ 


mi 


m2 -m3 J 




Iimil2m2 



{-lyi-h+rn,/ 21:^ + 1 ( ^1 ^2 '3 ^ . (3.18) 

^ \^ TOi 7712 — m3 J 



The Wigner's 3j (and, consequently, the Clebsch-Gordan) coefficients are real-valued, they are differ- 
ent from zero only if mi + ?Ti2 + rn^ = and h < Ij + Ik for all i,j, fc = 1, 2, 3 {triangle conditions), and 
they satisfy the symmetry conditions 

; /lyi+h+h ( ^1 ^2 ^3 

I —7711 —1^2 —7713 



h 


h 


^3 


mi 


7772 


"73 


h 


h 


^3 


mi 


7772 


7773 



I 7772 "^3 '77i 

where tt is a permutation of {1,2,3}, and sign{TT) denotes the sign of vr. ft follows also that for 7771 = 
7772 = "m-s — 0, the coefficients Cj^q; q are different from zero only when the sum li + 12 + I3 is even. Later 
in the paper, we shall also need the so-called Wigner's 6j coefficients, which are defined by 

a b e \ ^ V^ / 1^e+/+e+0 ( a b e\ f c d e \ f a d f \ f c b f 
c d f j- 2^^ '> \ape\-i5-e\a5-<p\^Pcj, 

(3.19) 
see |31| . chapter 9 for analytic expressions and a full set of properties; we simply recall here that the 
Wigner's 6j coefficients can themselves be given an important interpretation in terms of group represen- 
tations, namely they relate different coupling schemes in the decomposition of tensor product into direct 
sum representations, see [7] for further details. 

For future reference, we also recall some further standard properties of Kronecker (tensor) products 
and direct sums of matrices; we have 

®r=i iA^B^.) = m=lA^) i®7=iB.) , (3.20) 

{®7=iAi)®B = ©^^i(^,(^B) (3.21) 

and, provided all matrix products are well-defined, 

{AB (g)C) = {A(g) In) {B®C) . (3.22) 

Here, ©"^^A^ is defined as the block diagonal matrix diag {Ai, ..., A„} if Ai is a set of square matrices of 
order Vi x ri, whereas it is defined as the stacked column vector of order (X]"=i ''O x 1 if the Ai are r^ x 1 
column vectors. 

4 Characterization of polyspectra 

4.1 Four general statements 

The following result is well-known. As it is crucial in our arguments to follow and we failed to locate any 
explicit reference, we shall provide a short proof for the sake of completeness. Note that, in the sequel, 
we use the symbol a;, to indicate the (2/ + l)-dimensional complex-valued random vector defined in (|2.4p . 



Lemma 4 Let T be a strongly isotropic field on S'^ , and let the harmonic coefficients {aim} be defined 
according to \2.'/3^) . Then, for every I > and every g G SO{3), we have 

D\g)ai.=ai. ,1^0,1,2,... . (4.23) 

The equality Ii4-23\ l must be understood in the sense of finite- dimensional distributions for sequences of 
random vectors, that is, {4-2^ takes place if, and only if, for every fc > 1 and every < h < I2 < ■■■ < Ik, 

{D'Hg)ai,.,...,D'>'{g)ai,.}^{ai,.,...,ai,.}. (4.24) 

Proof. We provide the proof of (|4.24p only when k — I and ^1 = ^ > 1. The general case is obtained 
analogously. By strong isotropy, we have that, for every I > 1, every g e SO (3) and every xi, ..., a;„ € S'^, 
the equality (|2.ip takes place. Now, (|2.ip can be rewritten as follows: 

Kim I m ) K I m I m ) 

= j XI 51 "'™ 51 D^m^m{9)Ylm'{xi), ..., ^ ^ a,„ ^ D^^, ^{g)Yi^, [Xn) \ 
K. I m m' I m m' ) 

= l'^'^ahn'Yl,n'{xi),...,'^'^aijn'Ylm'{Xn)> , (4-25) 

L / m' I m' ) 

where we write 

alr,^' ^^al„^D'^,„^{g), (4.26) 

and we have used 

{Yi.m{gxi),...,Yi„^{gXr,)} = \Y.^'n.'mi9)Yi,rAxi), ...,J2Dl,Jg)Yi„,,{x,,) I . (4.27) 

l^ T7i' m' ) 

which follows from the group representation property and the identity p.9p . To conclude, just observe 
that (|4.25p implies that 



aim' = / T{gx)Yira' {x)dx, m' = -l,...,l, 

yielding that, due to strong isotropy and with obvious notation, Sj. — ai,. The conclusion follows from 
the fact that, thanks to (j4.26p . 

ai, = D^{g)ai.. 

■ 

The next theorem connects the invariance properties of the vectors {a;.} to the representations of 
50(3). We need first to establish some notation. For every < /i, ^2, •■•, ^n, we shall write 

Ai,...,„ ^ / {D'^[g)®D'-{g)®...®D'-[g)]dg , (4.28) 

JS0{3) 

^h...iA9) = D''[g)®D'-{g)<i^...®D'-{g), <? G 50 (3) , (4.29) 

and use the symbol Si^,,,i^ (whenever is well-defined), as given in formula (|2.5p . We stress that A/j^ ;^ 
and A/j ;^ {g) are square matrices with {2li + 1) x ... x (2/„ + 1) rows and 5;j...;„ is a column vector with 
{2li + 1) X ... X (2/„ + 1) elements. The following result applies to an arbitrary n > 2: see [TS] for some 
related results in the case n — 3, 4. 
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Proposition 5 Let T be a strongly isotropic field with moments of order n > 2. Then, for every < 
li,l2, ■■■,ln and every fixed g* G SO (3) 

Ai,...i„Si,...i,^ = Si,...i„ (4.30) 

Ai,..,A9*)Si,..,„ = Si,..,„. (4.31) 

On the other hand, fix n > 2 and assume that T{x) is a not necessarily isotropic random field on the 
sphere s.t. sup^ (£'|r(a;)| ) < oo. Then T {.) is P-almost surely Lebesgue square integrable and the nth 
order spectral moments of T exist and are finite. If moreover \4.S0^ holds for every < h < ■ ■ ■ < In, 
then one has that, for every g € SO (3), 

E [L>'i {g)ai,. ®---® D^"{g)ai^] = E [ai,. ®---® aij , (4.32) 

and T is n-weakly isotropic. 

Proof. By strong isotropy and Lemmas one has 

E {D'^{g)ai,,^ ...^ D^^^{g)ai^,} ^ E {ai,. ^ ... (E) ai,^} for all .g G 50(3) , ?i, ...,;„ G N". 

Now assume that g is sampled randomly (and independently of the {a;.}) according to some probability 
measure, say Pq, on 50(3). From the property p.22p of tensor products and trivial manipulations, we 
obtain (with obvious notation and by independence) 

E {D^'{-)ai,, (g) ... ® D^-{-)ai„,} = E { [D^'{-) (g) ... ® £»'"(•)] [ai,. ® ... (g) ai„.]} 

= Eo {D^'i-) (g ... (g £»'"(•)} -B{ah. ^ •■• «> a;„.} • 

Now, if one chooses Pq to be equal to the Haar (uniform) measure on SO (3), one has that 



En{D''{-)<i^...<»D'-{-)}^Ai,, 



■In 



thus giving (|4.30p . On the other hand, if one chooses Pq to be equal to the Dirac mass at some g* £ SO (3), 
one has that 

Eo{D''{-)(g>...(ED'-i-)}^Ai„„i,^{g*), 

which shows that (|4.3ip is satisfied. 

Now let T satisfy the assumptions of the second part of the statement for some n > 2. We recall first that 
the representation (|2.2p continues to hold, in a pathwise sense. To see that the nth order joint moments 
of the harmonic coefficients aim are finite it is enough to use Jensen's inequality, along with a standard 
version of the Fubini theorem, to obtain that 



E\aijn\ — E 



T{x)Yim{x)dx 



<E \Tix)r\Yimix)rdx 



< \ snp \Y,mixT \\ sup E\T{x)r 
2/ + l^"/' 



It is then straightforward that, if 5;^..j„ satisfies (|4.30l) . one also has that for any fixed g £ SO{3) 

E { [D'l (g) (g ... (g D^"{g)] [ai,. (g ... (g a,„.]} 
= [D'i(5)(g...(gi:)'"(g)] E[ai,. g)...(ga/„.] 
= [D''{g)<g>...<g>D'-ig)]Ai,,„i„Si,,„i„ 

= \[D''{g)®...®D'"{g)] f {D'^ (g) ® ... ® D'"ig)} dg\ Si,...i„ 

[ JS0{3) J 

{D'^igg) <» D'^igg) ^ ... <» D'-{gg)} dg\ Si,,„i^ 

50(3) J 

Azi..j„5/,...i„ = E{ai^, g) ...g)a/„.} , 
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which proves the n-th spectral moment is invariant to rotations. The fact that T is n-weakly isotropic is 
a consequence of the spectral representation 



Note that relation (|4.30p can be rephrased by saying that, for a strongly isotropic field, the joint 
moment vector E {ai^, ® ai^. (E) ... (E) ai„.} niust be an eigenvector of the matrix (|4.28p for every n > 2 and 
every < Zi < • • • < Z„. A similar characterization holds for cumulants polyspectra. Recall the notation 
Sf I introduced in (|2.6p . 

Proposition 6 Let T be a strongly isotropic field with moments of order n > 2. Then, for every < 
li,l2, ■■■,ln and every fixed g* G SO (3), 

^h...KSf^...i„ = Sf^„,^ (4.33) 

A/,...i„(g*)5r^..,„ = 5r,..,„. (4.34) 

On the other hand, fix n > 2 and assume that T{x) is a not necessarily isotropic random field on the 
sphere s.t. sup^, (-E|T(a;)| ) < cxd. Then T {.) is P-almost surely Lehesgue square integrahle and the nth 
order spectral moments of T exist and are finite. If moreover {^.33^ holds for every < h < ■ ■ ■ < In, 
then one has that, for every g £ SO (3), relation Il4-32\ l holds, and T is n-weakly isotropic. 

Proof. For every xi, ...,a;„ € S"^, write Cum{T(a;i) , ...,T(a;„)} the joint cumulant of the random 
variables T (xi) , ..., T (a;„). By using isotropy, one has that, for every g e SO (3), 

Cum{r(a::i) ,...,r(a;„)} == Cum{r(ga;i) , ...,T (gxn)} . (4.35) 

Hence, by using the well-known multilinearity properties of cumulants, one deduces that (with obvious 
notation) 

Cum{T(a::i) ,...,T(a;„)} 

= X! CvLin{ai^rni,-;ai^rnAYhmi{xi) ■ ■ -Yi^rn^iXn) 

lirrii ,...,1^171^ 

= X! Cum{ai^rni,-;ai^rnAYhmi{gXl) ■ ■ -Yi^rn^igXn) , (4.36) 

li7ni .....Ijinin 

and relations (|4.33p - (l4.34p are deduced by rewriting (j4.36p by means of the identity 

{Yl,rnA9Xl),■■■,Yl^„^A9Xn)} ^ lY.D'^,^^ig)Yi,rn'ixi),...,Y,Dt'^Jgm 

\ m' m' ) 

The second part of the statement is proved by arguments analogous to the ones used in the proof of 
Proposition m ■ 

We now present an alternative (and more involved) characterization of the cumulant polyspectra 
associated with an isotropic field. Given n > 2 and a partition tt = {hi, ...,hk\ G n(n), we build a 
permutation v"^ = {v'^ (1) , ...,1;'^ (n)) E S„ as follows: (i) write the partition 

TT = {5i, ...,bk} = {(*}, ...,*ffc,|) , ..., (it -.^ffc,!)} (4.37) 

(where \bj\ > 1 stands for the size of bj) by means of the convention outlined in Section [5] (that is, order 
the blocks and the elements within each block according to the lexicographic order); (ii) define v'^ = 6„ 
by simply removing the brackets in (|4.37p . that is, set 

v^ = (v^ {!),..., v^{n)) - (ij,...,iffc^|,z2,...,if,^|,...,zj,...,if,^l) . 
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For instance, if a partition tt of {1, ..., 6} is composed of the blocks {1, 3} , {6, 4} and {2, 5}, one first writes 
■K in the form tt = {{1, 3} , {2, 5} , {4,6}}, and then defines w'^ = (t;'^ (1) , ...,1;'^ (6)) = (1,3,2,5,4,6). 
Given n > 2, < ^i < • • • < ^„, and tt e 11 (n), we define the matrix 

^l..U = / {^'"'''' (5) ® ^'^"<^' (5) ®-® ^'^"'"' (ff)} rfff , (4.38) 

"'SO(3) 

obtained from the matrix A/j^.j^ in (|4.28p . by permuting the indexes U according to w^. Plainly, if v'" is 
equal to the identity permutation, then /S.J ^ = A/j...;,^. We also set, for every fixed g E SO (3), 

K...i„ (5) = ^'""<^' (5) ® ^'""<^' (5) ^ - ® ^'"'<'" (5)- 

Proposition 7 Let T be a strongly isotropic field with finite moments of order n > 2. ForO<li,l2,...,ln, 
define Sf_^ ;^ according to i2.6\) . Then, for every < h, h, ■••, ^n, a?irf every g £ SO (3) 

Sf,.,..,„ = E (-l)'"'(fc-l)!Ar^..,^i?[®,e6iaiJ^---^i=^[^^efc.ai.] (4.39) 

7r={bi,...,6fc}en(n) 

51 (-l)''"'(fc-l)!Ari...i„(<?)^[®»efcia/..]<8)---(8)S[®.e6,a/..] (4.40) 

7r={bi,...,6fc}en(n) 

On the other hand, fix n > 2 and assume that T{x) is a (not necessarily isotropic) random field on the 
sphere s.t. sup^. {E \T(x)\"') < c». Then, the nth order spectral moments and cumulants of T exist and 
are finite. If moreover l4-4(^ holds for every < ^1,^2, ■•■,^n and every g G SO (3), then one has that T 
is n-weakly isotropic. 

Proof. Fix TT — {bi,...,bk} G Il{n). By strong isotropy and Lemma SI one has that, for a fixed 
g* e SO (3), the quantity 

E [(S>ieb^D^' (5) ai^] (g>---(g>E [0.e6,D'* (g) a/,.] 
= Al_i^ (g*) E [®,g6ia,..] (D ■ ■ ■ (^ E [®iefc,a,,.] 

does not depend on g* , so that 

E [^if^b^ai^] ®---®E [®i(zb^ai.] 
= Af^ ;^ (5*) E [®,e6ia/..] ®---®E [^iefefca;..] 

= / Af^ l^{g)E[(E)ieb^al^]®■■■®E[®,(,b^au]dg 

JsO{3) 

= A^^i^E [^reb^a^.] (E) ■ ■ ■ (E) E [®jgb,a;^.] . 

To prove the second part of the statement, suppose that T{x) verifies sup^ {E |r(x)|") < 00, and that its 
associated harmonic coefficients verify (|4.40p . Then, for every fixed rotation g* S 5*0(3), 

7r={&i,...,6fc}6n(n) 

^ (_lf-i(fc_l)!x 

7r={6i,...,6fc}en(n) 

X [D'-'d) (g*) (g) • • • (g) d'-'m (g*)]E [(g>^eb^al^.] (g) ■ ■ ■ (g) E [(g^ebfcfli,.] 

7r={bi,...,bk}en{n) 

Yl {~l)'"^^ik-l)lE[(g>,eb^ai^,](g)---(g>E[(g>,eb^ai^,] 

7r={6i,...,6fc}en(n) 
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By the definition of cuniulants, this last equahty gives that 

E[D''^{g*)ai,. «) • • • «) D'"(.g*)a/„.] :^E[ai^,(^-- • (g) a/„ 
Since g* is arbitrary, the n-weak isotropy foUows from (|2.2p . ■ 



Remark. By combining (|4.33|) and ()4.39p we obtain for instance that the nth cumulant polyspectrum 
of an isotropic field verifies the identity 

Y. (-l)''"'(^-l)!Ar,...i„^[<»»e6ia/,.]®---®S[®,g6,a/,.]. 

5 Angular polyspectra and the structure of /^i^...i^ 

5.1 Spectra of strongly isotropic fields 

Our aim in this section is to investigate more deeply the structure of the matrix Ajj...;^ appearing in 
(|4.28p . in order to derive an explicit characterization for the angular polyspectra. As a preliminary 
example, we deal with the case n = 2. 



Proposition 8 For integers li,l2 > 0, one has that 



^hh 



SO(3) 



{D^Hg) ® ^'^5)} dg = 6\:iCl%iCl%y 



(5.41) 



that is: if li ^ I2, then Aj^/^ is a {2li + 1) (2/2 + 1) x (2/i + 1) (2^2 + 1) zero matrix; if li — I2, then 
Ai,i,^Ai,i, tsgivenbyC°XXC?°,J- 

Proof. Using the equivalence of the two representations D^^ (g) D^^ (g) and 0^/ _; |-C'^((?), as well 



as the definition of the Clcbsch-Gordan matrices, we obtain that 

/ {D'Hg)®D'^ig)}dg = Ci,i., f K=+|',U|^'(5)| ^5 

JsO(3) JsO(3) ^ ' ' J 



'x=\h-h\-' 



Q 



iib 



(5.42) 



Now, if ^1 7^ ^2, then the RHS of (|5.42p is equal to the zero matrix since, as a consequence of the Peter- 
Weyl theorem and for A 7^ 0, the entries of -D'^(-) are orthogonal to the constants. If h = I2, then the 
integrated matrix on the RHS of (j?^ becomes JgQ,^^ {^^xlo^^id)} ^9, that is, a (2Zi + 1)^ x (2/i + 1)^ 
matrix which is zero everywhere, except for the entry in the top-left corner, which is equal to one (since 
Iso(3) dg ^ 1)- The proof is concluded by checking that 



C/1/1 



SO(3) 



[<LoD\9)}dg 






Remark. Recall that Cf"; is a column vector of dimension (2/i + 1) (2^2 + 1), corresponding to the 
first column of the matrix Ci-^i^. Also, according e.g. to [31j formula 8.5.1.1], one has that 



^00 _ ] ( ^)\ l2X-m2 



mi= — ii,...,ii;m2 = — ^2 



Proposition [8] provides a characterization of the spectrum of a strongly isotropic field. 
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Corollary 9 Let T be a strongly isotropic field with second moments, and let the vectors of the harmonic 
coefficients {a/.} be defined according to 112. 2\) . Then, for any integers h^h > 0, one has that 

E{ai,.^ai,.} = |hll!^5;2^-™2c,^| (5.43) 

for some Ci-^ > depending uniquely on li. 
Proof. According to (|4.30p . one has that 

implying that E{ai^, (g) ai^,} is (a) equal to the zero vector for Zi ^ I2, and (b) of the form Ci^ x Cf^^^ , for 
some constant Ci-^, when Zi = 12- To see that C;^ cannot be negative, just observe that ai-^o is real- valued 
for every li > 0, so that (|5.43p yields that 

Ci, = {2h + l)xE (alo) ■ 

■ 

In the subsequent two subsections, we shall obtain, for every n > 3, a characterization of A;^...^;,^ and 
E{ai^,^ • • • (g) a/„.}, respectively analogous to (|5.4ip and (|5.43p . 

5.2 The structure of A^^ /^ 

We first need to establish some further notation. 

Definition B. Fix n > 3. For integers li, ...,Z„ > 0, we define Ci^,,,i^ to be the unitary matrix, of 
dimension 

n n 

n (2z, + 1) X n {21, + 1) , 

connecting the following two equivalent representations of SO (3) 

D^^{.)®D^^{.)®---®D^^{.) (5.44) 



and 



5Ai = |b-Zi| ®\2 = \h-\i\- ®a"„-i = |z'„-A„_2| ^ " '(■)■ (^-^^^ 



Remarks. (1) Fix /i, ..., Z„ > 0, as well as 5 G SO (3). Then, the matrix 

®tX~i.\ <hL-a,- <:'1T;„-a„_.| D'-'i9) (5.46) 

is a block-diagonal matrix, obtained as follows, (a) Consider vectors of integers (Ai, ..., A„_i) satisfying 
the relations I/2 - h\ < M < h + h, and |/fc+i - Afe_i| < Afe < h+i + Afe-i, for k = 2, ...,n - 1. (b) 
Introduce a (total) order -(q on the collection of these vectors by saying that 

(Ai,...,A„_i)-(o(A'i,...,A:,_i), (5.47) 

whenever either Ai < X'l, or there exists k — 2, ..., n — 2 such that Xj = A' for every j — 1, ..., fc, and 
Afe+i < X'k^i- (c) Associate to each vector (Ai, ..., A„_i) the matrix Z?-**"-! {g). (d) Construct a block- 
diagonal matrix by disposing the matrices D^^^^ (g) from the top-left corner to the bottom-right corner, 
in increasing order with respect to -<o. As an example, consider the case where n = 3 and li = I2 = I3 = 1. 
Here, the vectors (Ai, A2) involved in the direct sum (|5.45p are (in increasing order with respect to ■<q) 

(0,1), (1,0), (1,1), (1,2), (2,1), (2,2) and (2,3), 
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and the matrix (|5.46p is therefore given by 

( D^{g) 

1 



\ 



DHg) 



D'ig) 



DHg) 



(5.48) 



D' (g) J 



where the dots indicate zero entries, and we have used the fact that D^ (g) = 1. 

(2) The fact that the representation ()5.45p has dimension YVj^i ("^h + 1) is a direct consequence of 
formula (l37T4l) . 

(3) The fact that the two representations (I5.44|) and (|5.45l) are equivalent can be proved by iteration. 
Indeed, by standard representation theory, on has that (|5.44p is equivalent to 

which is in turn equivalent to 

By iterating the same procedure until all tensor products have disappeared (that is, by successively 
replacing the tensor product D^>'{.) ® 0^"+^ {■) with ^''^^^''^J'n^ 2-aJ^^'(-) ^°^ ^ = 2,...,n- 1), one 
obtains the desired conclusion. 



For every n > 3 and every Zi, ..., /„ > 0, the elements of the matrix Ci^,,,i^, introduced in Definition 
B, can be written in the form C, ^Z "i il^"^^. The indices (TOi,...,m„) are such that ~h < nij < L 
{i = 1, ...,7i) and label rows; on the other hand, the indices (Ai...A„_i,/i„_i) label columns, and verify 



the relations 1^2 — 'i| < Ai < /i + Z2, |^/c+i — A 



fc-i 



— 1) and — A„_i < 
A*n-i ^ A„_i. It is well known (see e.g. [21j) that the quantity C; ^^ Tm!""^ *^^^ ^^ represented as a 
convolution of the Clebsch-Gordan coefficients introduced in Section 



< Afc < Ik+i + Afe_i (fc = 2, . 

I 

namely: 



C, 



Ai,...,An,_i,/l„_i 
ll7ni ...InTUji 






E 



.c 



^n-2l-l-n-2 



^71-2 I Ml---A'7i-3 



E^Xlfll /-^A2/^2 y^A^_2M77-2 ^ 

Iimil2m2 \l^ilh■m3^^^^Xn-3^tn-3l■n-lmn-l ■ 



c 



A„-iM„-i 

A„_2M7i-2'nn7n 



A„-lM7l-l 

A„-2Mn-2'77 777„ 



Remark. Given an enumeration of the coefficients C^ ^^' "; m'""^' ^^^ matrix Ci^,,,i^ can be built 
(analogously to the case of the Clebsch-Gordan matrices of Section [3?^ by disposing rows (from top to 
bottom) and columns (from left to right) increasingly according to two separate total orders. The order ^^ 
on the symbols {mi, ..., to„) is obtained by setting that (ttii, ..., r7T,„) ^^ ("^'ij •■•j fri-'n) whenever either mi < 
m[, or there exists k = 2, ..., n—1 such that rrij = m', for every j — 1, ..., fc, and m,k+i < mi'j^,^. The order 
-<c on the symbols (Ai...A„_i, /i„„i) is obtained by setting that (Ai...A„_i,(U,„_i) -<c {^'l■■■^'n-l^^^'n-l) 
whenever either (Ai, ..., A„_i) ^o (A'l, ..., ^'n-ij j as defined in (|5.47p . or A^ = A^ for every i = 1, ..., ti — 1 
and /x„_i < /ij^_i. 

One has also the following (useful) alternative representation of generalized Clebsch-Gordan matrices. 
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Proposition 10 For every n > 3 and every li, ..., Z„ > 0, one can represent the matrix Ci-^,,,i^, as follows 



Q 



'h...U = {Chl2h...lr.-1 f^hu + l} [{®'C=\h-h 



J„+A„_3 



-/i|---^A„_2 = |Z„-A„_3| 



c 



A„-2i„ 



where Im indicates a m x m identity matrix. Also, one has that 



Ci^...i„ = [Ci.i^^hi 



3 + 1 



('p'2+il 



■>b+h 



-'Xi = \h-h\ 

where x stands for the usual product between matrices 



/2/„ + l)x (etl^/.I^A^a) 



^2;„+i 



i,.+A„_3 

A„_2 = |i„-A„_3| ^^^^-2lr^ 



Definition C. For every n > 3 and every li, ...,ln > 0, we define Ei^,,j^ to be the H"^]^ {21 j + 1) x 

117=1 (2^j + 1) square matrix 



E, 



h...i„ 



j2+il 

^Ai = |b-/i|- 



/„ + A„_2 rO T 

5a„_i = |/„-A„_2| '^A„_i^2A„_i + 1- 



(5.49) 



In other words, Ei-^,,,i^ is the diagonal matrix buih from the matrix (|5.46p . by replacing every block of 
the type D^"-^ {g), with A„_i > 0, with a (2A„_i + 1) x (2A„_i + 1) zero matrix, and by letting the 1x1 
blocks D° {g) — 1 unchanged. For instance, by setting n = 3 and h = h = h = ^ (and by using ([57 



one obtains a 27 x 27 matrix Em whose entries are all zero, except for the fourth element (starting from 
the top-left corner) of the main diagonal. 

The following result states that the matrix A/j...;^ can be diagonalized in terms of Ci-^,,,i^ and Ei-^,,,i^. 
Proposition 11 The matrix A;j..i^ can he diagonalized as 

(5.50) 
where Ei-^,,,i^ is the matrix introduced in Definition C. 

Proof. One has that 

D'l (g) (E) D''' {g) (g) ■ ■ ■ (g) D^" (g) dg 



^ii...in — Ci^,,,i^Ei^,,,i^Ci^^ 



Hx...l„ 



30(3) 



(5.51) 



SO(3) 



Ci^...i^ 



12+11 

^Ai = |;2-/i| '^A2 = |i3-Ai|- 



J,.+A„_: 



-A„_2i ^'""n5)Q:...„ 



dg. 



By linearity and by the definition of the integral of a matrix-valued function, one has that the last line 
of (|5.5ip equals 



Ci 



h...i„ 



^Ai=|;2-/i| ^A2=|;3-Ai| 



^A„_i=|;,.-A„_2| 



/ D^--ig)dg 

JS0(3) 



c, 



ll...l„ 



Now observe that, if A„_i > 0, then JgQ/^) E)^"^-^ {g)dg equals a (2A„_i + 1) x (2A„_i + 1) zero matrix, 
whereas Jqq(o\ D^{g)dg — Jgnr^) lcf.9 = 1. The conclusion is obtained by resorting to the definition of 
Eii...i^ given in (|5.49p . ■ 
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5.3 Existence and characterization of reduced polyspectra of arbitrary orders 

Combining the previous Proposition with ([5]), we obtain the main result of this paper. 

Theorem 12 If a random field is strongly isotropic with finite moments of order n > S, then for every 
li, ...,ln there exists two arrays -P;i....i„(Ai, ..., A„^3) and P^ i„(-^ii ■■■! A„_3), with \l2 — li\ < Ai < ^2 +'i, 
1^3 - All < A2 < ^3 + Ai, ..., \ln-2 - A„_4| < A„_3 < ln-2 + A„_4, such that 

Eai,ra....ai^^„ = i-ir- J2 ■■■T.^tn[^XtZ:.Ph....iA^i,---,K.-3) (5.52) 

X1 — I2 — I1 A„_3 

Cnm{ai,^„...,ai„.^J = (-1)"" £' •■• E ^w:.7„'^iZ":i^U....z.(Ai, ..., A„_3) (5.53) 

Ai— ^2— ^1 A„_3 

^Ai...A„_3;i„-m„ _ \ ^ \ ^ ^Ai/^i ^A2P2 pin.-mTi /c rA\ 

iimi....i„_im„_i / y •" / y Iimil2m2 Ximlsnis'" A„_3^„_3i„_im„_i ' \ ' / 

Remark. For afixed ri > 2, the real-valued arrays {Pjj...;^ (•): ?i, ..., ^„ > 0} and {P;*^ ; {■) '■ h, ■■■,ln > 0} 
are, respectively, the reduced polyspectrum of order n — 1 and the reduced cumulant polyspectrum of 
order n — 1 associated with the underlying strongly isotropic random field. 

Proof of Theorem 1121 We shall prove only (|5.52p . since the proof of (|5.53p is entirely analogous. 
By Proposition [5] and Proposition [Til if th^ random field is isotropic, then 

Jh...U — ^h■■■l,^-t^h■■■l,^^ll...l„^h■■■U : 

that is, because Ci^,,,i^ is unitary 

It follows that Si-^,,,1^ is a solution if and only if the column vector Q* ; Si-^,,,i^ has zeroes corresponding 
to the zeroes of Ei-^,,,i^, whereas the elements corresponding to unity can be arbitrary. In view of the 
orthonormality properties of C;* ; , this condition is met if, and only if, Si-^,,,i^ is a linear combination 
of the columns in the matrix Cj* ; corresponding to non-zero elements of the diagonal Ei^,,,i^. These 
linear combinations can be written explicitly as 

I2-I1 l3 + \l i„+A„_2 

I2-I1 i3+Al i„+A„_2 ( 1 

— Z^ Z^ •■■ Z^ \ 2^ ^hmil2m2^\itiil3m3---^>^n~2^tr,^2-lr,m„"l f ^ 

\l=l2-li \2=h-^l >^n-l=ln->^n-2 (_A'l---Mn-2 J 

X -Pii....i„(Ai, ..., A„_3, A„_2) 
h-h h+M '"+^" 

\ ^ J \ ^ /^Ai/^1 /^A2M 

ll3m3---^\n-2lin,-2-ln,m„ 



/ y ■'■ / y I / y Iimil2m2 Xil-iihms'"^^ 



Ai— ^2 — h A2— ^3 — Ai A„_i— /„ — A„_2 lMl---Mn-2 

X P(i....i„(Ai, ..., A„_3, A„_2)- 

Recalling again that 

f~iOm __ \ I X^2f— rK2X0 

^Iimil2m2 9/4-1 '1 ™i ™ ' 
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(see [3T], 8.5.1.1), we obtain that 



Ai— ^2— U A2— /s — Ai A„_i— /„— Ati_2 Imi---M7i-2 



EY^ Y^ J Y^ (-rAipi ^A2M2 I '-J dn A-m- 

/ ^ ■" / ^ I / ^ Iimil2m2 Ai^iiams ■■■ 2^ i ^ A„_2 ^11-2 

X Pli....l„{^l, ■■■, A„_3, A„_2) 



i2-ii i3+Ai /,7+A„_2 I j 

_ Y^ Y^ \^ < \^ /-fAi^ii ^A2M2 r^ln-rrin I _\\">-n \ y- 

/ ^ / ^ '" / ^ I / ^ iimii2m2 Ai/iii3m3 ■••^A„_3/j„_3.i„_im„_i V / J 

Ai=i2 — ii A2=i3 — Ai A„_i=i„ — A„_2 ^_ A'i---A'ji-2 j 

X P/,...j„(Ai,...,A„_3) 

I2-I1 i3+Al /„ + A„_2 

= 2^ Z^ •■■ Z^ C';^'„;; "7^^_"J""^P/i....i„(Ai,...,A„_3), 

Ai— /2~^1 A2^/3 — Ai A„_i^/ti~ Ati_2 

where we have set 

-Pii....i„(Al, ..., A„_3) := —— — F;j....;„(Ai, ..., A„_3,^n) ■ 

All there is left to show is that the coefficients of this linear combination are necessarily real. To see 
this, it is sufficient to specialize the previous discussion to the case where ra\=rfi2 — ■■■ — rUn — 0,, and 
to observe that, in this case 

.EaiiO---ai„o = 2_/--- zZ ^hoLiZ-iO -Pii---i,>(Ai, •■•, A„_3) 

Ai A„_3 

is real by definition (note indeed that the columns of Ci-^,,,i^^ are linearly independent). ■ 
Let us illustrate the previous results by some more examples. 
Examples. For n — 3, Theorem [T2] implies that, under isotropy 

^"'Iimi"'l2m2"'hm:i — \ >- ) '~'l-^mil2m2^^^2^3 ■ 

From this last relation, we can recover the so-called reduced bispectrum, noted bi-^i^i^, defined for instance 
in [15] , [20] and [21] , which satisfies indeed the relationship 



p _, ^/30 / (2^i + l)(2^2 + l) 
^hhh - 0Hl2l,^l^0l^0\^ (2^3 + 1)4^ ■ 

For n = 4 (i.e. the trispectrum, [15] ) we obtain the expression 

h+h 
i^O'limxO-l2m20'hm3ai 4,7714, = (^1) 2_^ ^lrmil2m2l3milil2l3liW 

\=\l2-li\ 
I2+I1 A 

^ Z_^ 2^ ^hmil2m2^\lJ.l3m3-^hl2huW ■ 

A=|i2-ii|M=-A 

The next result gives a further probabilistic characterization of the reduced bispectrum. 

Proposition 13 Fix n > 2. A real-valued array {Aij...;^ (•) : /i, ...,Z„ > 0} is the reduced polyspectrum 
of order n—1 (resp. the reduced cumulant polyspectrum of order n—1) of some strongly isotropic random 
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field if, and only if, there exists a sequence {Xi : / > 0} of zero-mean real-valued random variables such 
that 

l>0 

and, for every /i, ..., Z„ > 

E (Xi, • • • X,J = J2 - E CiX'C!o'Ah--ui>^i^ -. A„-3) (5.55) 

X1—I2 — I1 Ati-3 



(resp. 

Cum{X;,,...,X(J= Yl ••■E^i''^V;.t:fo'%....i.(Ai,...,A„_3).) (5.56) 



h+l 



2+11 



^-,Ai...A„_3/„0 

Ai— ^2 — ^1 A„_3 



Proof. We shall only prove (|5.55p . For the necessity it is enough to take Xi = aio, where aio is the 
harmonic coefRcient of index (/, 0) associated with a strongly isotropic field with moments of all orders. 
For the sufficiency, we consider first the (anisotropic) random field 

Zix)=J2XiYio{x). 

l>0 

Then, by taking T (x) — Z (gx), where g is sampled randomly with the uniform Haar measure on 5*0 (3), 
one obtains a random field with the desired characteristics. ■ 

There are two very important issues that are left open by Theorem 1121 As a first issue, it seems 
natural to look for characterizations of the reduced polyspectra P/i...i„, at least under natural models 
of physical interest. As a second point, we note that the explicit expressions provided in Theorem 1121 
depend on the ordering /i, ..., Z„ we chose for the decomposition of A/j...;^. In the next two sections, we 
try to address these (and other) points. 

6 Some Explicit Examples 

In this section we provide explicit computations for the reduced polyspectra Pi-^...i^ {n > 2), or Pf^ ^ , 
for some models of physical interest. Of course, the Gaussian isotropic fields can be easily dealt with. 
Indeed, in this case one has that Pf ^ — for all n > 3. In what follows, we shall therefore be concerned 
with polyspectra of Gaussian subordinated isotropic fields, that is, random fields that can be written as a 
deterministic and non-linear function of some collection of Gaussian isotropic fields. In general, this class 
of random fields allow for a clear-cut mathematical treatment, whilst covering a great array of empirically 
relevant circumstances. 

6.1 A simple physical model 

The general Gaussian-subordinated model has the form 



T = Y1 fi^i [TgI^E (T2) = /iTg + /2(t2/S (t2) - 1) + ..., (6.57) 

where fj is a real constant, Hj{.) denotes the jth Hermite polynomial (see e.g. (SHDj and Tq is a 
Gaussian, zero-mean isotropic random field. Note that we have implicitly defined the sequence of Hermite 
polynomials in such a way that Hi [x) = x, H2 {x) = x^ — 1, H3 (x) = x'^ — 3x, and so on. In this section. 
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when no further specification is needed, the spectral decomposition of the underlying Gaussian field Tq 
is written 

Ira 

We shall sometimes use the following notation 

q 
T = ^ai„rYi„-,^^fjai„i{j)Yijn, (6.58) 

Im J — 1 



aimU) = / ^ Hj ITg (x) I^E (Tg) 1 Yi^ {x) dx, (6.59) 

aim = ^ aim (j) ■ (6.60) 

i=i 

For instance, models of Cosmic Microwave Background radiation are currently dominated by assump- 
tions such as the Sachs- Wolfe model with the so-called Bardeen's potential (see e.g. [5] or [H]). The 
latter can be written down explicitly as 

T^Ta + fNLiT^-ET^) , (6.61) 

where /jv^ is a nonlinearity parameters which depends upon physical constants in the associated "slow- 
roll" inflationary model (see e.g. [5]). Note that (|6.6ip has can be written in the form (|6.57p . by setting 
/i = 1, /2 = /nl X E {Tq) and fj = 0, for j > 3. The value of the constant fNL x E (T^) is expected to 
be very small, namely of the order 10"'^ [5]. To simplify the discussion, we now assume that ETq — 1. 
In this case, by using (|6.58p - (|6.60p . one has that 

aim — aim + fNLaim{2) , 



(2) = / T^Yimdx = / V V ae,m^ 



a;m(2) = / T Yimdx = / /^/^ at^miai^m2Ylim.iYi^m2Y imdx 



C2 

"^ lif,2 mim.2 



_ V- V- „ „ / (2^1 + 1) (2^2 + 1) ^,0 r'm 

— Z^ Z^ "^i™i"«2m2y ('9/-l-lU7r '-'^lO£20^<?imif2m2 



(21 + 1Utt 

ili2 mim2 I ^ ' 

It follows that 

Ci := E\aim\'^ = Ci + 2/^^ 2_, C'iiC/2 4j^(2l + I) (^'"o'2o) ' 

'l'2 



SO that 



l^an^) - Z^^^^'-Z^^^^' + 2/wL2^^h^i2 7^^T2 Z^(^ii0/20J 

I I hh ^ ' I 

1 

= Var{TG)^ !ljyaT{H-2ijG)) , 



V2^ + V+9f2 \S-C (^il±l)l 



as expected, due to the orthogonality properties of Hermite polynomials. For the bispectrum, we obtain 
therefore 

Eai^miai2m2ai3m3 = E{{ai^mi + /2a;inii(2))(a/2m2 + /2a;2m2(2))(az3m3 + /2a/3m3(2))} 

= f2Eai-^^mi\^)ai2m2ai3m3 + f2Eai-^miai2m2('^)ai3m3 

+f2Eai^miai2m2ai3m3 (2) + /2 Eai^mi (2)a/2m2 (2)a;3,„3 (2) 

_ (_T\m3(-il3-m3 p 
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where 



Phhh = 6/^ J ^^^'(2^3 +^1)4^ ^^ ^''o° 20 {ChCi, + C/,C,3 + Ci^Ci,} (6.62) 

, .3 V r'lo r'^0 r'30 (2^1 + i)(2^2 + i)(2^3 + 1) ^^ ._„. 

^i^2^3 V(4'r)' 



V2/3 + 1 I '3 '2 ^1 

The lack of symmetry with respect to the ^3 term is only apparent and can be easily dispensed with by 
permuting the multipoles in C^i'^i „ or using expression (|3.18p . Formula (|6.62p is consistent with the 
cosmological literature, where (|6.63p is considered a higher order term and hence neglected (see again 

([IS]))- 

6.2 The Connection with Higher Order Moments 

We now provide a simple result, connecting the reduced polyspectrum with the higher order moments of 
the associated spherical random field. 

Proposition 14 The following identity holds for every isotropic field with finite moments of order p and 
with a reduced polyspectrum {Pii...ip (■) : h, ■■■, Ip > O}; for every x G S'^ , 



{2li + 1) ■ ■ ■ {2lp + 1) \^ p, 1^ \ ^^Ai...Ap_3ipO 

i„-20 ■ 



l\...lp V Ai...Ap_3 



Proof. We use the trivial fact that 

I2I + I 



T{x) = T(0) = Y. ««>lo(0) = Y. °'o 



47r 
; I 



where is the North Pole and we used the fact that, for to 7^ 0, Yim (0) = and Yiq (0) — \/^j^ (see 
e.g. [23 Chapter 5]). Hence, 



'(2/i + l)---(2Zp + l) 



li...lp V 



Y^ / (2/i + 1) ■ ■ ■ (2/p + l) v^ p /^ ^ ^^Al...Ap_3^pO 

li ...Ip V Ai ...Ap_3 



Example. Take T = HgiTc), where iJ, is the gth Hermite polynomial. Then ETP = Cp, {E'T^}'^^^^, 
where Cpq € N denotes the number of Gaussian diagrams without flat edges with p rows and q columns 
(see [30]). Therefore, one has the identity 



E {2h + l)...(2^p + l) sr p (\ \ ^^Al...Ap_3^pO 

li...lp V Ai...Ap_3 
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6.3 The xl polyspectrum 

Previously in (|6.63p . we have implicitly derived the "xf bispectrum" , that is, the bispectrum associated 
with a field of the type T = H2 (Tq), where Tq is Gaussian, centered, isotropic and with unit variance. 
More precisely, with the notation (|6.58P " (|6.60p . one deduces from (|6.63l) that 



Eai.rn, {2)ai,m, (2)a;3„,3 (2) (6.64) 

2_^ WiO£20^flMl^2A'2^f30£40'-'£3A'3^4M4^f50f60'-'f5M5^ 



^1^2^3X Ml---A'6 

/ (24 + 1)(2£2 + 1) (2^3 + 1)(24 + 1) (24 + 1)(24 + 1) ^, 
^Y (2/i + l)47r (2/2 + l)47r (2/3 + l)47r 

X E {a^ i^ii 0^2/^2 a^3M3 ^<?4M4 ^^SMS 0^6/^6 } 

- sr_l V3-m3 V- Wio ^i.o ^Z30 (2^1 + 1)(24 + 1)(2€3 + 1) 






7(4^ 

^SSUa ^: \:]{C.C,M. (6.65) 

see f3P, p. 260 ; p. 454]. We now wish to extend these results to polyspectra of order p = 4,5,6 for 
random fields of the type T — H2{Tq), where (as above) Tq is Gaussian, centered, isotropic and with 
unit variance . As anticipated, here we focus on cumulants instead of moments. We have the following 
result. 

Proposition 15 The cumulant x (diimi (2) , ..., aipm^ (2)) (p = 4, 5, Q) associated with the harmonic co- 
efficients of an isotropic random field of the type H2 (Tg) (where Tq is Gaussian and isotropic, with 
angular power spectrum {C; : / > 0} j given by 

X (ai^m, (2) , ..., a,^„,^ (2)) = (-1)'.—. ^ C^n^::^^^. x ^'^ (^i' -' ^p-^) ' 

Ai ...Ap_3 

where the reduced cumulant polyspectrum < Pf^ ; (•) : li, ..., ^p > > is given by 



p'^'i f\\ — A9i I (2A + 1) v^ /^ r' r^'io r^'ao r^'*" r^'^o 
\ii+i2+i2+ii J ^1 '2 A 1 J A h k 



x(2£, + l)...(24 + l)(-l)'-'— M I I I \\ ,; ,; ,; \ forp^A 



pC;l 1^ \ \ - Q«/| / (2Ai + 1) (2A2 + 1) V^ ^ ^ ^liO W2O ^(40 W5O ^(30 „ 
Ui...h\'^i^'^'2) - -^S^y (47r)5(2;5 + 1) Z^ '-'«i---'-'''5'-'«i0£20'-'«20£30'-'«30£40'-'£40£50'-'«50£i0X 



pC:i /x x X \ _ on^n / (2Ai + 1) (2A2 + 1) (2A3 + 1) >r-^ i^Q i^Q jgO /-./sO r'^sO r<i-i.o 

■f/i...i6l^l''^2,A3j -^»4Uy (47r)6(2Z + l) ^ '-"'^i---'-'^6Wi0^20'-"^20<?30'-"<?30£40'-'(?40<?50'-'<?50£60'-"£60<?i0X 

X(2., + 1M24H-1)(-1)— -{^'^ I ^'KJ' I ,^^^}{,^; I U'^--^"'"- 
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Proof. The result can be proved by means of the standard graphical techniques for convolutions of 
Clebsch-Gordan coefficients, as described in [211 Chapters 11 and 12]. Here, we only provide the complete 
proof for the case p — 6. Let {aim} be the random harmonic coefficients associated with the underlying 
Gaussian field Tq. By definition, the field H2 (Tg) admits the expansion 



;>o m=~l 



where 



aim {2) = 22 ^il-miai2m2 Yi^mi {x) Yi2m2 {x) Yim {x)dx 

iimii2m2 

aiimiai2m2 X (-1) X 

Y mi 7712 ^TTl I 

iimil2m2 






h 


^2 



r 


I \ /(2£i + l)(2^2 + l)(2/ + l) 





JV 47r 




Im /(2^1 + 1)(2£2 + 1) 



47r (2/ + 1) 



limi(.2in2 

By using once again the multilinearity of cumulants, one obtains that 
Cum {ai^mi (2) , ■■■,ai^ms (2)} 

= / ^ ■■■ / ^ Cum |a^jj„jiia^i2mi2 7 •■•J ^feimei'^^62m62l X 

<?iimii<?i2mi2 ieimeileimei 

6 f 

11 ) ^<'jim3l^j2mj2'-'^3l0£320 



,,.„.,. /(2£,i + l) (2^,2 + 1) 



47r (2;^ + 1) 



Foragivenlm = (^1177711,^1277712; ...;4i"76i,42?n62), the quantity Cum {a^j^m^^a^j^m^^, ..., a£5imeia£e2>n62} 
is computed as follows: 



• Build the 6x2 matrix 



A(lm) 



^1177711 ^1277712 

^2lW21 i22m22 

^3l"73i 42TO32 

£4l"74i ^2^42 

4l77l5i 42W52 

4l"761 42?T762 



• Define the class M (A (Im)) of connected, Gaussian non-flat diagrams over A, that is, every 7 £ 
M (A (Im)) is a partition of the entries of A (Im), into pairs belonging to different rows; moreover, 
such a partition has to be connected, in the sense that 7 cannot be divided into two separate 
diagrams. For instance, an element of M (A (Im)) is 

7 = {{^1177711, £21»772l} {^22^22, 42*7732} {41^31, 4l»774i} 

{42TO42, 42»7752} {41^61, 4l»776l} {42TO62, ^12»77l2}} 

• For every 7 G M (A (Im)), write 



(where 6^ is the usual Kronecker symbol) 



24 



• Use the standard diagram formula (see again [301 )j to obtain that 

Gnm{ae^^,niiaei2mi2i ■■■i^i6imei^i62m62} — /^ '^ (7) • 

7GM(A(lm)) 

It follows that 

Cum{a;^„j (2),...,a;g™g (2)} 

Im 7eA/(A(lm)) J = l K V \ 1 J ) 

where the first sum runs over all vectors of the type Im = {inmn, £i2mi2; ■•■;^6i"T'6ij^62'7i62)- The proof 
now follows directly from graphical techniques. In particular, the previous term can be associated with 
an hexagon, having in each vertex an outward line corresponding to a "free" (i.e. not summed up) index 
liirii, i = 1, ...,6. An expression for convolutions of Clebsch-Gordan coefficients corresponding to such a 
configuration can be found in [3H p. 461], eq. 12.1.6.30. From this, standard combinatorial arguments 
and a convenient relabelling of the indexes, we obtain that 



{n'=i (2A, + 1)} 
p-,(A„A„A3) = 384o^Vm:TTr^^-'^ 



Xi+es+le+h 



/ . (2^1 + 1) • • • (2i() + l)Cfj...6fg6^^Q£^Q6^^Q^^Q6£^Q^^Q6^;^Q£^Q6£^g^gQ6^gQ£^Q 

l...«6 

h h h ]\ U h h ]\ 4 4 h ]\ 4 4 h 

Note that 3840 — 2p~^ {p—l)l = 2^5! is the number of automorphisms between graphs belonging to 

M(A(lm)). ■ 

We recall that the Clebsch-Gordan coefficients {C'^q^q} ^^^ identically zero unless a + 6 + c is even; it 
is hence easy to see that the previous polyspectra are non-zero only if the sum {^1 + ... + Ip} is even as 
well. 

From the previous Proposition, we can derive the corresponding expressions for the cumulant polyspec- 
tra for xi random field. 

Definition B. We say the random field T^2 has a chi-square law with v > \ degrees of freedom if 
there exist v independent and identically distributed Gaussian random fields Ti such that 

„ law rr,2 . , rr2 

^Xl = ^1 + ■■■ + -l-u ■ 

It is trivial to show that T-^2 is mean-square continuous and isotropic if Ti is. We have the following 
Proposition 16 The cumulant polyspectra of T^2 (for p > 2) are given by 

Proof. Note that the cumulant polyspectra of order p > 2 of T~^2 coincide with those of the centered field 
r^2 — ET^2 (due to the translation- invariance properties of cumulants). Then, the proof is an immediate 
consequence of Proposition 1151 and the of the standard multinearity properties of cumulants. ■ 
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7 Further Issues and Applications 

The purpose of this final Section is to introduce what we view as promising directions for further research, 
where the ideas of this paper may perhaps yield further insights. We shall delay to future work a more 
thorough investigation of the issues which are left open below. 

7.1 Representations of the Symmetric Group 

As a further link between representation theory and higher order angular power spectra, we mention 
the following. It is to be stressed that the decomposition of A/^.j^ that we achieved in the previous 
Proposition 111! is by no means unique. In particular, what we did was to choose a particular sequence 
of "couplings" , i.e. we partitioned tensor products of the Wigner's matrices _D' in a specific order before 
decomposing them into direct sums. Alternative partitions yield different eigenvectors and therefore, 
different expressions for the polyspectra/joint moments . Alternatively, we could maintain the same 
coupling scheme (for instance, "start always from the first pair on the left" , as we did earlier) but acting 
on (li, ...Tin) by the symmetric group S'„. However, not all coupling schemes can be achieved by simply 
permuting the elements of {hjh, •■•,^n)- This is the well-known problem of parentheses in Mathematical 
Physics (see for instance [7]). 

We suggest here that one can establish a link between alternate expressions for the angular polyspectra 
and representations of the symmetric group. More precisely the alternate expressions that we find for 
the polyspectra Pij....;^^(Ai, ..., A„_3) of a strongly isotropic field (with n-moments) must be such that, 
for every permutation tt G &„, 

Z^--- Z^ ^iimi..".7,!irmr-"l^'i----'"('^l'---''^"-3) 
Ai A„_3 

■■■ Z^ '-'7r(ii)mi....Tr(i„_i)m„_i^'^('i)----'r(i,i)l^li •■■'^„-3J ■ 

Now let us multiply both sides by C; ^^" ";"^^ "^ "_ , where (A", ..., X'^_^) is fixed, and sum over (toi, ...m„). 
In view of the unitary properties of Clebsch-Gordan coefficients we obtain for the left-hand side 



Ai A„ 

on the right-hand side we get 



E^Ai...A„_3;(„-m„ I ^ y^ „Ai ... A„_3;i„-m„ „ ,. , 

mi...m„ I Ai A„_3 

EY^ J Y^ ^A'/...A^_3;i„-m„ Ai...A„_3;i,>-m,> p /\ \ \ l 

■■• Z^ ] Z^ ^lirni....l„-imn-i'~'lrmi....l„-im^-i-^h....lr.\M,--;^n-3) > 

Ai A„_3 Lmi...m„ ) 



E^Ai...A„_3:(„m„ I V^ V^ ^Ai...A„_3:(„m„ , , , - 
mi...m„ ^ A'j A^_3 

EY^ Y^ „A'/...A^_3;i„m„ A'i...A;_3;i„m„ , , , , (7 C7\ 

■■■ Z^ Z^ ^7r((i)mi....7r((„_i)m„_i'-^7r(ii)mi....7r(i„_i)m„_i ^('i)----'r('")^^l' ■■■' ^ri-3^ I'-"'; 



Similarly as in the previous section, the sum of products of Clebsch-Gordan coefficients on the right hand 
side can be expressed in terms of higher order Wigner's coefficients. Since this section is just informal, 
for brevity's sake we do not give explicit expressions (see e.g. ^TJ Chapter 10]). The two expressions 
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(j7.66p and (j7.67p imply that, for every fixed (Zi, ...,/«) and every permutation tt, there exists a square 
matrix A ((/i, ..., Z„) ; tt) such that 

Ph....l„ = A{{li, ...,ln) ;7r}P7r(ii)....7r(i„), 

where Pi^...i^ is the vector with entries Pi^...i^ (Ai, ..., A„). We conjecture that in this way one can build 
a representation of the symmetric group (3„ on the vector space generated by admissible polyspectra 
Pii....i„- If this is indeed the case, some important questions are left open: for instance, whether or not 
the representation is faithful (see |10j). and whether these ideas can lead to algorithms for the numerical 
simulation of representation matrices, along the lines of what we shall pursue in the next subsection. 

7.2 Random data compression 

In this subsection we shall show how we can exploit the previous results to develop a probabilistic 
algorithm to compress information on Clebsch-Gordan coefficients. Note first that 

# {C&™. ■■ hMM < L, \cl%Z.r...\ ^ 0} « OiL^) ; 

it is therefore clear how for most applications the storage of Clebsch-Gordan coefficients for future usage is 
simply unfeasible, whatever the supercomputing facilities (for instance, for CMB data analysis, L w 3x 10'^ 
is currently required, so that the number of Clebsch-Gordan coefficients to be saved would exceed 10^°). 
Let us consider again a chi-square field as defined before, i.e. 

T^2(x) = H2{Tg{x)) = ^a/„(2)l,„(x) ; 
we have proved earlier in (|6.65p that 

where 

h -s V^ rho W20 ^M (2£i + l)(2^2 + l)(2^3 + l) 1 f £1 £2 ^3 \}r r r X 

which can be calculated analytically and stored, with storage dimension 

*[hiM, : hMM < L, \C\%^^\ ^ 0} « 0{L^) . 

Let us assume we simulate B times r^2(x), which is trivially done by simply squaring a Gaussian field: 
the latter is obtained by sampling independent complex Gaussian variables with variance C;. We store 
the triangular arrays {a;„}, , ^ . __, ; ' * ~ f ' •■•' -^5 "tieve the dimension is of order B x L^ . We can 
then recover any value C,^!"^, ^ by means of the Monte Carlo estimate 



Iimil2m2 hhh / ' 



(i) (i) (i) 

"' - '^I2m2"'hm3 



B 

which requires B steps and B x L^ + L'^ storage capacity, as opposed to L^ storage capacity by the direct 
method. We leave for further research a more thorough investigation on the convergence properties of 
this algorithm; we stress, however, that the procedure we advocate is completely general, i.e. it does not 
depend on peculiar features of the group 5*0(3) we are currently considering. We believe, hence, that 
similar ideas can be implemented for the numerical estimation of Clebsch-Gordan coefficients for other 
compact groups of interest for theoretical physicists. We leave this and the previous issues in this Section 
as topics for further research. 
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